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In this work, we generalize the paraxial ray-tracing formulas to include nonparaxial rays. For a refracting (reflecting) 
spherical surface, a new single meridional formula is derived. This formula can be easily reduced to a paraxial formula. 
It can also be applied to any aspheric (or general) surface with a known equation. Also, a new exact ray-tracing 
procedure for a centered system of spherical surfaces is derived. In this procedure, we apply just two simple equations 
for each surface of the system, which, to the best of our knowledge, makes it the shortest analytical ray-tracing 
technique ever. This procedure can be applied in some other applications. For example, it can be reduced to a 
new single paraxial formula that can be easily used to trace a paraxial ray propagating through a system of spherical 
surfaces. Also, it is applied to derive an exact meridional formula for both thick and thin lenses that can also be 
reduced to a new paraxial formula different from the Gaussian one. These results led us to easily derive an exact 
formula for the longitudinal spherical aberration for both thick and thin lenses and also for a single refracting 
(reflecting) spherical surface. Numerical examples are provided and discussed. © 2017 Optical Society of America 
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1. INTRODUCTION 


In the literature of geometrical optics, there are several exact ray- 
tracing procedures suggested [1-7]. In most cases, any one of 
these procedures is presented as a series of equations that is ap- 
plied to each surface of an optical system. For many decades, these 
procedures were continuously developed and widely used in the 
field of optical design. At the present time, there are several high- 
quality, ray-tracing computer programs available to the optical 
designer. The ultimate goal of optical design is to reduce ray aber- 
rations to get the sharpest image possible [2,3]. However, practical 
applications should not be the only task of exact ray tracing, 
Another possible task is to help us understand the behavior of 
light rays in the nonparaxial region. In Gaussian optics, we have 
a group of simple linear equations describing the behavior of 
paraxial rays in optical systems, but we really miss a similar group 
of exact equations for nonparaxial rays. All optical engineers 
and optics students know these paraxial formulas and use them 
frequently, but they never see their generalized forms. Each of 
these generalized forms should be reducible to the paraxial form 
in a few simple steps. Mathematically, it also should be con- 
tainable in the body of any theoretical work concerning optical 
design or any other topics. The present work aims to derive 
some of these exact nonparaxial formulas using simple mathe- 
matics. We think that deriving such formulas can be achieved 
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successfully if we avoid using the concept of principal planes 
in the paraxial region. We replace the principal planes by the cur- 
vature centers of the optical surfaces. For example, in the case of 
refraction at a spherical surface, we determine the positions of the 
object and the image by measuring the distances between each of 
them and the curvature center of the surface. For the case of a 
thick lens, the locations of the object and the image are deter- 
mined by measuring the distances between each of them and 
the curvature centers of the first and second surfaces of the lens, 
respectively, and not to the principal planes. This way, the ob- 
tained meridional formulas can be easily reduced to simple para- 
xial forms. In the present work, our attention is mainly focused 
on meridional rays and spherical surfaces. The aim of the work is 
not to design a ray-tracing procedure for computer software and 
commercial purposes, but to explore how such general formulas 
can be derived, how they will seem, and how they may help as an 
introduction to a theory for nonparaxial optics. 


2. SIGN CONVENTION 


(1) Light rays travel from left to right. 

(2) The radius of curvature, R, is positive for a convergent 
(refracting/reflecting) surface and vice versa. 

(3) The distance S, between the object and the vertex V of 
the surface is positive when the object is located to the left of V 
and vice versa. 
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(4) The distance Sı between the image / and the vertex V 
of the surface is positive when the image is located to the right 
(left) of the refracting (reflecting) surface and vice versa. 

(5) The distance K, between the object and the center of 
curvature C for a (reflecting/refracting) surface is positive if the 
object is to the left of C and vice versa. 

(6) The distance K, between the image and the center of 
curvature C for a (reflecting/refracting) surface is positive if the 
image is to the right of C and vice versa. 

(7) The angle that the ray makes with the optical axis is 
positive if the latter has to be rotated in the counterclockwise 
direction (through the acute angle) to coincide with the ray and 
vice versa. 


3. REFRACTION AT A REFRACTING SURFACE 
A. Description 


Figure 1 shows a spherical surface separating two media of re- 
fractive indices 7, and 74. Let us assume that O is a point object 
and J is its image formed by the surface. Let K, and K, be the 
distances of both O and J measured from the center of curva- 
ture C. In the present work, all distances are measured to the 
centers of curvature of the refracting (reflecting) surfaces. This 
greatly simplifies the mathematical work. The incident ray OP 
makes an angle 0, with the optical axis, while the refracted ray 
PI makes an angle 0). 

Our target now is to calculate the location of the image 7 
(i.e., to calculate the distance K,). In geometrical optics, we are 
used to solving this problem by using a series of trigonometric 
equations. In this work, we will prove that the distance K} can 
be calculated using the formula 


1 
ra L sin sin'a- sin! 9, , (1) 
where 
K, sin 8 
a decd 2 
a = (2) 


Equation (1) is a new exact meridional formula for refraction at 
a spherical surface. It can be easily used to calculate the image 
distance K, by direct substitution, as will be shown in a 
numerical example. The derivation of Eq. (1) is based on 
the geometry of Fig. 1. Let g and @ be the angles of incidence 
and refraction, respectively. According to the sign convention 
mentioned above, the angle 0, is positive, while 0, is negative. 
So, we can write 


Fig. 1. Meridional refraction at a spherical surface. 


a=0,-0,=9-8 (3) 
where æ is the deviation angle. Therefore, 
0,=90,-9+8. (4) 


For simplicity, we can take n = n/n, where n is the relative 
refractive index and nı > n,. So, we can write sin y = n sin 8. 
Now let 


a = sin Q. (5) 
Then, g = sin! a and 0 = sin"! a/n. Thus Eq. (4) can be 


written as 
0, = 0, - sinta + sinta/n. (6) 
In Fig. 1, we have 


. ri K,sin@, 
sin Q R R (7) 
Therefore, 
K, sin 0 
ae 8 
aa Kes 8) 
And from the triangle CP” we have 
Kı r2 


= ——, 9 
sin 90 sin(-0,) (9) 
The relation between r; and 7, can be easily found as follows: 
sin g = rı/R and sin 0 = r2/R. In this case, sin p/ sin 0 = 
rı/r2 and by Snell’s law, this directly leads to ry = nr. 
Therefore, 


K, sin 8 
p= al cane (10) 
n 
From Eq. (9), we get 
1 n 
— = ——— sin(-0,). 11 
Kane (11) 


Finally, we substitute from Eqs. (6) and (8) into Eq. (11) to 
get Eq. (1). 


B. Numerical Examples 


(1) In Fig. 1, let us calculate the distance K, of the image 7 
from the center of curvature of a refractive spherical surface if 
K, = 18 cm, R = 8 cm, n = 2 and 0, = 15.825489. Using 
Eq. (2), we get a = 0.613593627. Then we substitute into 
Eq. (1) to easily get Kı = 33.849960 cm. 

(2) For the case of a concave (divergent) surface, let us 
consider the following numerical example: R = -8 cm, 
K,= 12cm, n= 2, and 0, = 8.783323. Substituting in 
Eqs. (1) and (2) we get a=-0.229047 and K = 
-3.439578 cm. According to the sign convention, K, must 
be negative, since the image is formed to the left of the curva- 
ture center. To calculate the distance S} of the image J from 
the vertex v of the surface, we just add R to Ky (i.e. 


Sı = Kı + R= -3.439578 - 8 = -11.439578 cm). 
Remarks: 
(1) Equation (1) is not only used to determine the position 
of the image, but it can also be reduced to the paraxial formula, 


as we shall see in the next section. It also can be used to get an 
exact spherical mirror meridional formula. 
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(2) This is right for both meridional and skew rays. For ex- 
ample, if the point O in Fig. 1 is located out of the optical axis 
and Op is a skew ray, then we can consider a secondary optical 
axis OC passing through any point C located at the normal to 
the surface. In this case, the plane OPC is considered as a sec- 
ondary meridional plane and Eq. (1) can be applied to get Ky 
of the image that will be formed at the secondary optical axis. 
The details of this idea will be published in a future work. 


C. Paraxial Approximation 


In the paraxial region, the ray’s inclination angles are very small. 


So, for any angle x we have sinx =x and sin™!x = x. 


Therefore, Eq. (2) reduces to 
a= —. (12) 


In this case, æ ~ 0. So, substituting into Eq. (1) we get 


1 n a 
— = ge 1 
Kı zje n a a 


1 n on-l 
Ki K, R’ 
This is the paraxial formula written in terms of the distances K, 
and K,. It also can be written in the form 


(14) 


1 4 n _n-l 
Sı-R S, +R R 
where the distances S, and S; are shown in Fig. 1. The two 


paraxial formulas (14) and (15) are equivalent to the 
Gaussian one, 


f (15) 


taioe 
S, S& R’ 


So, we have three formulas for the single spherical surface. 


(16) 


D. Parallel Incidence 


When the incident rays are parallel to the optical axis, then 
ri = K, sin 0, =A, (17) 


where / is the height of the point of incidence P from the op- 

tical axis. In this case, we find that a = A/R and Eq. (1) will 
have the form 

h h 

K, = i sin sin i sin“! a (18) 

This is the parallel incidence equation in the trigonometric 

form. The nontrigonometric form can be obtained by expand- 

ing the sin function in Eq. (18) and this leads to the formula 


E +R 
VRR -P - SR? 
In Eq. (19), the positive sign is chosen for the convergent 
refractive surface and vice versa. 


(19) 


Ky 


4. REFLECTION AT A SPHERICAL MIRROR 
A. Description 


In Fig. 2, O is a point object located at a distance K, from the 
curvature center C of a spherical mirror, and J is the image 


Fig. 2. Meridional reflection at a spherical mirror. 


formed by the mirror for this object at a distance K, from 


C. Since 
1 1 


K = sin @, (20) 
where 7 is given by r = K, sin 8, therefore, 
1 
K = RAT sin[0, + 2a]. (21) 
Now, 
sin a = (K, sin 0,)/R. (22) 
Thus, Eq. (21) can be written as 


1o 1 
K, K, sin @, 
This is the meridional formula of a spherical mirror for any 


angle @,. It can be used to derive another nontrigonometric 
formula for the mirror having the form 


1 1 2 h\2 
abe Sy 1 | =), 24 
K, Kı R 9 (ea) 
where / is the height of the point of incidence from the optical 
surface. To derive Eq. (24), we use (21) and (22) to get 


(23) 


K, sin 0 
sin fo, + 2sin"! Enej 


Ki = Riva <info, + 2a] B 
= z4 bin 4, cos(2a) + cos 8, sin(2a)] J” 


Using the identities 


cos(2a) = 1-2 sin? a 


sin(2a@) = 2 sin a cos a, (26) 
we get 
1 1 : : nig 
— = ——|sin 9, - 2 sin 9, sin* a 
K, Rsina 
+ 2 cos 0, sin a cos a], (27) 
1 sin 0, 2 sina 


K, ETTI 7 sin 0, sin a - cos 0, cos a]. (28) 


Using Eq. (21) we get 
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1 1 2 
— = — += cos . 2 
K, KTR cos[0, + a] (29) 


In Fig. 2, we have 


cos[O, + a] = cos y = y R? -/7/R. (30) 


By substituting into Eq. (29), we finally get Eq. (24), which is 
the simplest meridional equation for a spherical mirror. 


B. Another Derivation to the Mirror Equation 


Equation (24) can be derived following simpler steps. In Fig. 2, 
we see that pọ = a + y, and 0, = y - a. So, we have 


sin y = sin a cos y + cos @ sin y, (31) 


sin 0, = sin y cos @- cos y sin a. (32) 
Subtracting Eq. (32) from Eq. (31), we get 
sin ọ - sin 0, = 2 sin a cos y. (33) 


It is clear that sin gp = r/K, sin 0, = r/K sin æa = r/R, 
and cos y = \/1-4°/R*. So, substituting into (33), we 
get Eq. (24). 

Remarks: 


(1) It should be noted that Eq. (23) can be, directly, ob- 
tained from Eq. (1). In this case, we must remember that 
Eq. (1) was derived for a convex refracting surface. So, accord- 
ing to the sign convention, Ris negative for a convex mirror, 
and n = -1. Thus, Eq. (1) reduces directly to Eq. (23). 

(2) Equations (1), (23), and (24) can be applied also for 
skew rays. For example, if the point O in Fig. 2 is located 
out of the optical axis and Op is a skew ray, then we can con- 
sider a secondary optical axis OC passing through any point C 
located at the normal to the surface. In this case, the plane OPC 
is considered as a secondary meridional plane and Eqs. (1), 
(23), or (24) can be applied to get K, of the image formed 
at the secondary optical axis OC. An extension of this work 
to the case of skew rays is to appear in the near future. 


C. Numerical Examples 


In Fig. 2, let us calculate K, for the image formed by the con- 
cave mirror if we know that R= +6 cm, K, = 7.2 cm, 
0, = 22.5. This problem can be easily solved by direct sub- 
stitution in Eq. (23), where we get Kı = 2.825833 cm. If 
the height / is given instead of 0, substitution into 
Eq. (24) leads to the same result. 


D. Paraxial Approximation 


In the paraxial region, we consider that 0, ~ 0 in Eq. (23). So, 
sin 0, ~ 0,, and the quantity between brackets also tends to 
zero. Thus, Eq. (23) reduces to 


1 1 2K 0 
= ra ga. 4 
zn eae a 
So, we get 
-1 1 2 
K, K = R (35) 


Equation (35) is the paraxial formula of reflection at a spherical 
mirror. It can be also written in the form 


ee eee 
RsS, R-S R 


where S, and S; are the distances of both the object O and the 
image J from the mirror. So, Eqs. (35) and (36) are equivalent 
to the famous Gaussian formula, 


1 1 2 
oe 37 
S, S R e) 
So, we have three forms for the paraxial equation of a spherical 
mirror. It should be noted that Eq. (24) of the mirror can also 
be reduced to the paraxial form by just putting 4 = 0. 


(36) 


5. REFRACTION (REFLECTION) AT ASPHERIC 
SURFACES 


Equations (1) and (24) can be applied for the case of aspheric 
refracting (reflecting) surfaces. For example, in Fig. 3, a 
meridional ray is refracted at a surface with equation f(x, y) = 
0 (in two dimensions). If the height / of the point of incidence 
is known, we can use Eq. (1) to calculate K} as follows: 
The three quantities K,, 0, and A can be easily calculated 
as follows: solving f(x, 4) = 0 for x, then we can get 0, = 
cot !(x,//). The angle y is given by y = tan”! 1/ f'(x, 4). The 
length A, which corresponds to R in Fig. 1, is normal to the 
tangent at the point of incidence. So, we have A = // sin y and 
xı = h cot y. Therefore, the distance K, is given by 
K, =x, +x. Now, substituting in Eq. (1) for K,, 0,, and 
A we get K,. For the case of reflection, we follow similar steps 


using Eqs. (23) and (24). 


6. RAY-TRACING PROCEDURE FOR A 
CENTERED SYSTEM OF SPHERICAL SURFACES 


A. Description 


In the literature, there are several procedures for meridional ray 
tracing. In an exact ray-tracing procedure, we usually use a 
series of trigonometric (or vector) equations for each one of 
the successive optical surfaces. In this section, we describe a 
new simple procedure in which we use just two equations 
for each one of the surfaces. It should be stressed here that 
the technique only works in the case of spherical surfaces. 
We believe this procedure will be the shortest and fastest 
analytical technique published in the literature. 


Fig. 3. Meridional ray refraction at aspheric surface with equation 


fy) = 0. 
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B. Procedure Derivation 


In Fig. 4, we have two spherical refracting surfaces. The object 
O is located at a distance K, from the center of curvature of 
the first surface. Assuming that g, and @; (not shown in figure) 
are the angles of incidence and refraction, respectively, the 
deviation angle of the first surface is 


ô = 91 - Pi- (38) 
Now, let 6, and 6, be the angles made by the incident and 
refracted rays, respectively, with the optical axis. In this case, we 
have 
= 20) (39) 
Thus, 
9, =0, -pı + pi- (40) 


—¢ : Sa 1 — jp] mar 
Let a, = sin gj; thus, pı = sin” a and g = sin! =, 


: K, sin 0 
Also, we have to note that a4; = sin g; = a Thus, for 
the first optical surface, we get the following two equations: 


K, sin 0, 
4 = TR ’ (41) 
1 
0, = 0, - sina, + sin! Teli, (42) 
n 


1 
Equations (41) and (42) are sufficient for the first surface. We 
do not need to calculate the image position formed by this sur- 
face, because this will be done only for the last surface of the 
optical system. For the second surface (see Fig. 4), similar steps 
are followed. We should remember that 8, is negative accord- 
ing to the sign convention given in Section 1. Now, the 
deviation angle 6, is given by 


ô, = 9, - 02 = P2 - P% (43) 
where p, and g, are the angles of incidence and refraction at 
the second surface, respectively. So, we can write 

0 = 01 - Gr + 9}. (44) 


= ee | Po as 1 ya 
Let a, = sin 2, so we have p = sin™' a and p3 = sin’ 2%, 


and Eq. (44) can be written as 
0, = 0; - sina, + sin! alec (45) 
ny 
Figure 4 shows that 


r3 _ r3 +d, sin 0; (46) 


a = sin P2 = 


R R 


Fig. 4. Meridional refraction through a system of two spherical 
surfaces. 


_ (n,/m)r, + d; sin 0, a P dı sin an (47) 
R, mR R, 


4a 


where r} = (”,/m,)r}, as was proved in Section 3.A, 
and r; = K, sin 0,. 


So, 
nK, sinô, dı sin 0, 
=< 2 4+ l 48 
az mR R (48) 
Now, since K, sin 0, = a,R,, we can finally write 
1 |7, . 
k= —a,R, + dı sın 0, . (49) 
Ry |m 


If we need to calculate the distance K, of the final image 
formed by the second surface, we just substitute in the equation 


r4 
Ky= ; 50 
2 sin(-0,) Bo) 
This can be rewritten as 
1 na 1 . 
Se 05. 51 
Ky (*) aR ve wn 


So, we finally use Eq. (51) to calculate the image distance K3, 
by just substituting for 0, and a) from Eqs. (45) and (49). 
Equation (51) can be easily derived starting from Eq. (50) 
as follows: 


1 1 
K = TH sin 0z, (52) 
n= n= (ntd sin 8) =" aR) (53) 


Substituting from Eq. (53) into Eq. (52), we get Eq. (51). 
For an optical system with several surfaces, two equations 
having the forms of Eqs. (45) and (49) are used for each surface, 
and finally, the image position of the last surface is calculated. 
Summary of the procedure: 
For the mth surface of a centered system of spherical surfaces 
we have the equations 


1 
4n => [ee Ry aI d m1 sin Oni} (54) 
Rn An-1 

PES | . 1 Žm-1 
On = Oy - Sin! dp + sin! Z ap (55) 

Am 

1 n 1 
= -| — in 6p 56 
K m (=) AmyRin a ( 
Remarks: 


(1) For each surface of the optical system, we just use the 
first two equations, (54) and (55). The third equation, (56), is 
only used for the last surface. 

(2) In this procedure, we consider the distances d „, between 
the successive curvature centers instead of the distances £, sepa- 
rating these surfaces. So, all the distances d,, for the optical 
system must be given before applying the procedure. The fol- 
lowing equation can be used to calculate the distance d„ in 
terms of the separation ¢,, between any two successive surfaces 
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ami = bml" (Rm-1 a Rn)» (57) 


(3) For the first surface (m = 1), there is no need to use 
Eq. (54), since a, = Eate is already known because K,, 
0, and R; are all given. 


C. Numerical Examples 


Assume that we have a system of four centered spherical sur- 
faces that form a final image 74 for a real object O located on 
the optical axis of the system at a distance K, from the curva- 
ture center C4 of the first surface. Now, we need to calculate 
the exact distance Ky of the final image 74 if we know the 
following: 


0, = 17.309724, K;=22 cm, R, = 10 cm, 
R, = -8 cm, R; = 12 cm, R4 = -10 cm, 

ti = 5 cm, t = 5 cm, t3 = 8 cm, n=l, 
ny = 1.2 m = 1, nz = 1.5, ng = 1. 


To solve this problem, we first need to calculate the distances 
d, between the curvature centers using Eq. (57), where we 
easily get dı =-13 cm, d,=25cm, d3=-14cm, 
d4 = 3 cm. Also, for the first surface, the constant a, is given 
by Eq. (41) as a, = 0.654581. Now, for each surface we just 
substitute in Eqs. (54) and (55) to get the following: 


first surface, a4) = 0.654581, 0; = 9.479589; 
second surface, a. = -0.414224, 0, = 4.143784; 
third surface, a; = 0.481920, 03 = -5.926743; 
fourth surface, a4 = -0.530095, 04 = -26.583586. 


The final image distance K4 is given by Eq. (56) as 
K4 = 17.768432 cm. 


This numerical example proves the simplicity of this ray-tracing 
technique, and that it is the shortest in geometrical optics. In 
the previous calculations, all values are restricted to only seven 
digits. The higher the number of digits, the more precise is the 
result. 


D. Another Procedure 


In this section, we get another procedure in which we calculate 
the image distance for each optical surface. The procedure 
consists of two equations for each surface of the optical system. 
In Eqs. (54)—(56), and for simplicity, consider the relative re- 
fractive indices: p4 = n/n, Hy = m/n, and so on. Now, 
from Eq. (56) we get 


1 
SS ane td ain Oy: (58) 


Am tRm-t | — 
By = 1- À (59) 
lniRa K m-i 
Also, we substitute for 0,,_; from Eqs. (58) into (55) to get 
R 
0, = sin! a sin™! a,, - sin”! Silai (60) 
Hm Hin m-1 


Substituting from Eqs. (60) into (56) we get 


1 By Ra 
UmK mai] 
(61) 
Therefore, we have a new ray-tracing procedure in which we 
use just the two equations, (59) and (61). 


a 
sin |sin™la, - sin"! + sin 
m Am m m 


E. Paraxial Formula for a System of Spherical 
Surfaces 


In the previous procedure, we noted that the values of a,,_1, 
Am-2> Am-3» etc., can be obtained by Eq. (59). Now let us 
assume that 


d 
D= l= 2 
So, Eq. (59) can have the form 
R 
gom ae Da Dare Di (63) 
A m-14¢m 


Now, in the paraxial region, 0, ~ 0 and a, ~ 0 in Eq. (41). 
Thus using Eq. (63), we find that a,,~0 for all values of 
m. Therefore, in the paraxial region, Eq. (61) can be written as 


1 Hm Am Animi 
= = 7 64 
K A,R [en Hm Y aA ( ) 


Combining Eqs. (59) and (62), we get 


lmi Ryn 
== Dy: 
Am Hn1Rm m-1 (65) 
Substituting for a,,_, from Eqs. (65) into (64) we get 
1 Hm 1 Rn, 
= 1- i 
Ka Rn | Hm 7 ian (66) 


From Eq. (62) we get Dm-1Km-1 = [K m-1 -4ml and 
Eq. (66) can be rewritten as 


1 H Mm - 1 
= = 67 
"eu =K mi Rn ( ) 


Since My, = nm/nm-1> Eq. (67) can have the final form 


Nn Nin Am ~ Ron 

K + iw OR : (68) 
Equation (68) is the paraxial formula that can be used to trace a 
paraxial ray through a system of spherical surfaces and to cal- 
culate the position of the final image formed by the system. 
This will be shown in a numerical example. It should be noted 
that for the first optical surface (m = 1), the distance K, in 
Eq. (68) must be negative, since it represents the distance of 
a virtual image formed by an imaginary surface (m = 0) pre- 
ceding the first surface. Since the image is located to the left of 
the curvature center, K, is negative according to our sign con- 
vention. Thus, for the first surface, Eq. (68) can be written as 


(69) 


m m 


Ng ny _ 4 -% 


Ki K ç R 


F. Numerical Example 


Consider a centered system of four spherical surfaces. Let us 
calculate the final paraxial image J, formed by the system 
for a point object O placed at a distance K, from the curvature 
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center of the first surface. The following are given: 
K,=22cm, R= 10cm, ky =-8cm, Rk; = 12 cm, 
Rg = -10 cm, ¢4} = 5 cm, t7 = 5 cm, 43 = 8 cm. a, = 1, 
ny = 1.2, m = 1, n; = 1.5, n4 = 1. To solve this problem, 
we also need to calculate the distances d,,, between the curva- 
ture centers using Eq. (57), where we get dy = -13 cm, 
dy = 25 cm, and d; = -14 cm. Now we apply Eq. (68) 
for each surface to get the position of the image formed by this 
surface. So, we easily get the following: 


first surface, Kı = -28.947368 cm; 
second surface, K, = -31.825173 cm; 
third surface, K = 65.487885 cm; 
fourth surface, K4 = 40.087735 cm. 


(1) This method is different from that used in Gaussian op- 
tics. For an optical system of 7 surfaces, the final image position 
relative to the vertex of the last surface is given by 


S m 7 K m7 Rir 
(2) We can easily calculate the exact value for the longitudinal 
spherical aberration of a system of n spherical surfaces using the 
previous procedure, i.e., Eq. (68), and the exact procedure, 


Eqs. (54)-(56). 


7. THICK LENS MERIDIONAL EQUATION 
A. Description 


In Fig. 5, a point object O is placed on the optical axis of a thick 
lens at a distance K, from the curvature center C, of the first 
surface. The final image J formed by the lens is located at a 
distance K, from the curvature center C, of the second surface. 
The meridional ray OP makes an angle 0, with the optical axis. 

In this section, we prove that the exact meridional equation 
for this lens has the form 


1 m . |. B . nB ., m (B-aR, 
r = in | sin = sin -sin > 
Ky n,B j à mR mR nı dı 


d 
B= jar + ari sin (o, - sin"! a, + si! a) (71) 
Lil 


K, sin 0 
TEEN 72 
a Ry (72) 
Equation (70) can be reduced to the Gaussian (paraxial) 
thick lens equation, as will be proved in this section. To derive 


Eq. (70), we apply Eqs. (54) and (55) as follows: for the first 


surface (m = 1) we have 


Fig. 5. Propagation of a meridional ray in a thick lens. 


K, sin 0, 
a, = R, 2 (73) 
0, = 0, - sin™! a, + sin"! EPN (74) 
ny 

For the second surface (m = 2) we have 

1 /n, g 
4, = (Zar + di sın o); (75) 

Ry \m 

0, = 0, = sin! 4) + sin! Pi i (76) 


na 


The image location K is given by Eq. (56) as 


1 na 1 
=- in 05. 77 
Ky (=) aR ee 


Substituting from Eq. (76) into Eq. (77) we get 


1 = at) 
Ky 7 nyah, 


From Eqs. (74) and (75) we have 


sin 4, - sina, + sin! ta, (78) 
n 


d 
a= Po jar + A in (o, - sin™! a, + sin! *a)] 
n ny 
(79) 


Assuming that 


d 
B= jar ge EY mn (o, - sin”! a, + sin! ta)|, (80) 
o ny 
Eq. (79) becomes 
n 
= —? B, 1 
a nR, (81) 


Also from Eqs. (74) and (80) we have 


B-a R 
0 = sin! i = ‘)| (82) 
1 1 


Substituting from Eqs. (81) and (82) into Eq. (78), we 
get Eq. (70). 

B. Numerical Example 

For the thick lens of Fig. 5, we have the following data: 
R, = 10 cm, R, = -8 cm, 0, = 17.309724°, K, = 22 cm, 
n =1, m = 1.2, m =1 and d = -13 cm. To calculate 
Ky, we substitute into Eqs. (70)-(72), and we get 
a, = 0.654581, B = 3.976547, and Ky = -55.031208 cm. 
C. Parallel Incidence 


If the object O is at infinity (i.e, K, =œ and 0, = 0), 
Eq. (72) reduces to a; = 4/R, and the lens is described by 
the following equation: 


1 m . |. B . B . ,n,(B-h 
— = sin | sin - sin - sin ; 
Ky n,B mR mR nı di 
(83) 
d h h 
B= h + ae sin (sie = - sin”! alk (84) 


N.B.: If the lens is placed in air (i.e, n, = 1, n) =” and 
m = 1), Eqs. (83) and (84) will have the simpler forms 
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1 1 B B B-h 
K, = z sin sin "R - sin”! R - sin“! l 5 (85) 
where 


h h 
= E (eee wedi ay 
B= |} + nd, sin (sin z 7 sin *)|. (86) 


ni 
So, Eq. (85) is an exact meridional thick (or thin) lens equation 
in case of parallel incidence. 
D. Paraxial Approximation 


In the paraxial region, the exact lens formula, Eq. (70), reduces 
to the following paraxial form: 


1 n m ë m\ 1 nı -m 
= - - ’ 87 
K, 7 n,K ,E (2 2) RE mR, en 


where 
= h - ( z 1) AA a nel (88) 
n 


This can be proved by making use of the fact that in the paraxial 
region, 0, œ 0. So, from Eqs. (71) and (72), we get a, ~ 0 and 
B~0. Therefore, Eq. (70) reduces to 


l am |nB nB m (B-aR (89) 
Ky 7 n,B mR, mR nı dı f 


This can be rewritten as 


il 1 m n ak, 
K, R, (2-1) j (=) h- B } er 


From Eq. (72), we get 0, ~ a,R,/K,. So, Eq. (71) reduces to 


= l di mdi 
B = ak, h - (= = 1) R, + nel (91) 


Now, let us assume that 


E= h E (z z 1) a nel (92) 


Ng 


Then, Eq. (91) can be written as 
B= a RE. (93) 
Substituting into Eq. (90), we get 


nE (m-n n 
en a-ge- o 


Substituting from Eq. (92) into Eq. (94) and rearranging terms, 
we finally get Eqs. (87) and (88). 

N.B.: If the lens is surrounded by air (i.e., 2, = m, = 1 and 
nı = n), then Eq. (87) becomes 


1 1 n-1\fı1 1 
mt ree ( n eral as 
where 
= dı nd, 
E= fio a (96) 


E. Power of a Thick Lens 


In the literature, the power P of a lens is simply the inverse of 
the focal length. Now, if we let K, ~ œ in Eqs. (87) and (88), 


the distance K represents the focal length of the lens measured 
from the rear focus to the curvature center of the second sur- 
face. So, it can be replaced by the symbol f.. Thus, we have 


1 1 1 n nı -m 
= 2 2 P 97 
f. (= ~) RE ( mR ) ( ) 


pei m 


If the lens is surrounded by air (i.e., n, = m, = 1 and 7; = n), 


then Eqs. (97) and (98) take the simple form 
1 n-1 1 1 
f: ( n ) (ae -z) a) 


= (n- 1)d, 
E, = h aa | (100) 


So, combining Eqs. (99) and (100), we get 


1 j=l 1 1 
z=( n (re g) a 


Equation (101) is a new form for the paraxial thick lens equa- 
tion. In this equation, the quantity A can be considered as the 
power of the lens. In this case, we have to define a new power 
unit. In the Gaussian formula, the power is 1/ f „, where f, is 
the focal length measured to the principal plane. In a previous 
work [7], we proved that the relation between f, and f, is 


Se= Effy (102) 


N.B.: Eq. (101) can be derived from Eq. (68) in much sim- 
pler steps. 


where 


where 


8. LONGITUDINAL SPHERICAL ABERRATION 
A. For a Thick (Thin) Len 


The longitudinal spherical aberration (Long. SA), of a thin lens 
is measured by the difference 


(Long. SA) = S, - S,. (103) 


In Eq. (103), S, is the distance of the image formed by the 
paraxial ray measured to the lens, while S, is the corresponding 
distance of the image formed by a marginal ray incident at a 
height 4. Equation (103) can be written as 


1 1 
(Long, SA) = S,S,(—-—) = S,5,L- (104) 
S S 


In fact, the Long. SA can be minimized if the quantity L, is 
minimized. So, it is useful to derive an exact formula for 
the quantity Z, of a thin or thick lens. In the present work, 


distances are measured to the center of curvature, so 


Eq. (104) will have the form 


(Long, SA) = (Ky)p(Ko)s an i ra 
n p 


= (Ky) (Ko) Le (1 05) 
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1:6: 
1 1 


(K2), and (K2), are the distances corresponding to S,, and S,, 
respectively. Now we write Eqs. (85) and (101) as 


1 1 jin 1B ,B-h 


— sin - sin” — sin” 


(Kə), B nR, R, nd 


1 n-1 1 1 
= = , 108 
(Ky), ( n ) (x -(n-1)d, w) a 


Substituting from Eqs. (107) and (108) into Eq. (106), we get 


f (107) 


1 
L= 3 sin sin A - sin™! — - sin | 


n-1 1 j ? 
-( n (reor g) 


h h 
— = PPS i L inl 
B= Ç + nd, sin (sin 2 A Ii (110) 


nka 1 


where 


To the best of our knowledge, Eq. (109) is new, and it may have 
a possible theoretical importance. 


Remarks: 


(1) All the equations concerning a thick lens, including 
Eq. (109), can be applied for the case of a thin lens by just 
neglecting the thickness ¢ in the equation d} = ż - (R, - R2). 

(2) In the previous equations, the lens is assumed to be 
surrounded by air, and 7 is its refractive index. For the case 
of three refractive indices 1,, nų and m, as in Fig. 5, 
Eq. (109) can be easily generalized using Eqs. (83), (84), 
(97), and (98). 


B. For a Refracting Spherical Surface 


An exact formula for L, of a refracting spherical surface can be 
obtained using Eqs. (14) and (18) and also Eqs. (14) and (19) 


as follows: 


n-1 


h h 
L= ” sin sin g “| -R (111) 


h 


n-1 


1 
L=} VRR -P -VR-P|- (112) 
Equations (111) and (112) are two different exact formulas for 
L, of a refracting spherical surface. 


C. For a Spherical Mirror 


Using Eqs. (24) and (35), we can easily get an exact formula for 
L, of a spherical mirror as follows: 


2 BN? 2 
nati an 


9. CONCLUSION 


In this work, our target was to derive general ray-tracing for- 
mulas that can be reduced to simple paraxial equations. This 
could be easily achieved because the positions of the object and 
the images were measured relative to the curvature centers of 
the surfaces, not to the principal planes. The obtained formulas 
are new and can be easily used by optics students as generalized 
forms of the Gaussian equations. The work is just an approach 
toward a much more perfect theory for nonparaxial optics that 
may appear in the future. 


REFERENCES 


1. W. J. Smith, “Optical computation,” in Modern Optical Engineering 
(McGraw-Hill, 1990), pp. 308-404. 

2. D. Malacara and Z. Malacara, Handbook of Optical Design (Marcel 
Dekker, 2004). 

3. R. Kingslake and R. B. Johnson, “Meridional ray tracing,” in Lens 
Design Fundamentals (Academic, 2010), pp. 25-45. 

4. F.L. Pedrotti and L. S. Pedrotti, “Geometrical optics,” in Introduction to 
Optics (Prentice-Hall, 1993), pp. 34-36. 

5. P. Mouroulis and J. Macdonald, “Rays and foundations of geometrical 
optics,” in Geometrical Optics and Optical Design (Oxford University, 
1997), pp. 11-13. 

6. W. J. Smith, “Image formation: geometrical and physical optics,” in 
Handbook of Optics, W. G. Driscoll and W. Vaughan, eds. 
(McGraw-Hill, 1978). 

7. H. A. Elagha, “Exact ray tracing formulas based on a nontrigonometric 
alternative to Snell’s law,” J. Opt. Soc. Am. 29, 2679-2687 (2012). 


